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The Large-Scale Orientations of Disk Galaxies 
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ABSTRACT 

We use a 380 /i _1 pc resolution hydrodynamic AMR simulation of a cosmic filament to 
investigate the orientations of a sample of ~ 100 well-resolved galactic disks spanning 
two orders of magnitude in both stellar and halo mass. We find: (i) At z = 0, there 
is an almost perfect alignment at a median angle of 18°, in the inner dark matter 
halo regions where the disks reside, between the spin vector of the gaseous and stellar 
galactic disks and that of their inner host haloes. The alignment between galaxy 
spin and spin of the entire host halo is however significantly weaker, ranging from a 
median of ~ 46° at z = 1 to ~ 50° at z = 0. (ii) The most massive galaxy disks 
have spins preferentially aligned so as to point along their host filaments, (iii) The 
spin of disks in lower-mass haloes shows, at redshifts above z ~ 0.5 and in regions 
of low environmental density, a clear signature of alignment with the intermediate 
principal axis of the large-scale tidal field. This behavior is consistent with predictions 
of linear tidal torque theory. This alignment decreases with increasing environmental 
density, and vanishes in the highest density regions. Non-linear effects in the high 
density environments are plausibly responsible for establishing this density- alignment 
correlation. We expect that our numerical results provide important insights for both 
understanding intrinsic alignment in weak lensing from the astrophysical perspective 
and formation and evolution processes of galactic disks in a cosmological context. 

Key words: cosmology: theory, large-scale structure of Universe - galaxies: forma- 
tion, evolution - methods: numerical 



1 INTRODUCTION 



Numerical simulations and analytic calculations have shown 
that the gravitational growth of tiny density perturbations 
in the ACDM cosmological model leads to a wealth of struc- 
tures over cosmic time. The spatial distribution of gravi- 
tationally bound structures follows a web-like pattern com- 
posed of dense clusters, filame nts, sheets and otherwise near- 
empty void regio ns (see e.g. IShandarin fc Zeldovichl [l989l : 
iBond et al.lll99rj y This "cosmic web" is the weakly non- 
linear manifestation of the large-scale tidal field and its 
formation can be readily understood in the first order La- 
grang ian perturbation theory of the Zel'dovich approxima- 
tion (|Zel'Dovichl ll97Cf ) . Recent large scale galaxy surveys 
have observationally con firmed the presence of the cosmic 
web in the Universe (e.g. IColless et al.ll200ll ; iTegmark et al.l 
|2004) . 

Since galaxies form and evolve together with the cos- 
mic web in which they are embedded, it is a question of 
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high relevance whether the properties of galaxies depend on 
scales beyond the most immediate vicinity and thus reflect 
the effect of large-scale tidal fields on the assembly of the 
galaxy. General ellipsoidal density perturbations collapse se- 
quentially along their three axes leading to matter accreting 
first onto a sheet, then collapsing to a filament before matter 
finally streams along the filaments into the densest regions. 
As a consequence of these processes, the large-scale mor- 
phology of the cosmic web might leave its imprint on the 
orientations of galaxies through the anisotropy of accretion 
and mergers. 

Recently, evidence has accumulated that galaxies are 
indeed not oriented randomly with respect to one another. 
On small scales, several studies report a radial alignment 
of satellite galaxies towards the centre of mass of the 
group or cluster in which they reside (e.g. lYang et al. I l2006l ; 
Faltcnbachcr et al. 2008). Such an alignment is pre- 

dicted by N-body simul a tions of dark matt e r halos (e.g. 
Bailin fc Stemmed 120051 : IPereira et all 120081 : iKnebe et all 



2008). It is commonly believed that t idal torques are re - 



sponsible for this radial alignment (e.g. IPereira et al]|2008h . 
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as tidal forces on these relatively small scales (< 1 Mpc) are 
sufficiently strong to affect halo shapes on short time-scales. 

More surprisingly, however, observational results indi- 
cate that the orientations of galaxies are correlated up to 
scales of order ~ 100 times larger t h an the vi rial radii of their 
host haloe s. [M andclbaum ct al. (200jSj), lOkumura et al.l 
l|2009h and iFaltenbacher et all lj2009h detect an alignment 
of galaxy shapes for the luminous red galaxies (LRGs) in 
the S PSS survey out to the largest probed scales (~ 80 
Mpc). iHirata et all (|2007T ) report a similar result for red- 
shifts 0.4 < z < 0.8. For lower luminosity galaxies, in par- 
ticular fo r L < L t , no alignm ent has been detected yet 
(see also ISlosar fe White! [200!) i although some results us- 
ing N-body simulations indicate spin alig nment around voids 
l|Brunino et al.ll2007l ; ICuesta et alj |2008l. Correlations of the 
spin of galaxies have been the subject of several analytical 
models (see e.g. ISch afcr 2009, for a review, and references 
therein) . 

Compared with more local effects, the relative weakness 
on galaxy scales of tidal forces arising from the 3> IMpc 
scales requires that such effects maintain spatial coherence 
over rather long time scales for them to be able to affect 
galaxy properties. Also it is not clear how they are coun- 
teracted by local non-linear effects in the most immediate 
vicinity of the galaxy: non-linear gravitational effects as well 
as ram pressure exerted on the disk gas could potentially 
lead to a reorientation of galactic angular momentum. Fur- 
thermore, the impact of hot in contrast to cold mode accre- 
tion separated b y a halo mass of ~ 4 x 10 11 h^My of gas 
onto galaxies (c f. Birnboim & Dckcl 2003; Keres et al.ll2005l ; 
iDekel fc Birnboim|[200r3 ; lOcvirk et aTll2008h on the orienta- 
tion of the resulting disks is not clear. 

While understanding the formation of galaxies in the 
context of their host large-scale structure is a key ques- 
tion for galaxy formation theory, a precise knowledge of 
correlated galaxy orientations is of high importance also 
for future high-precision weak lensing studies of dark en- 
ergy. Weak lensing signals can be severely affected by 
systematic contaminations due to large-scale alignments 
in the spatial distribution of cosmic matter on scales 
from a few to more than a hun dred Mpc (d u e to g alaxy 
ellipticity-ellipticity alignment: ICatelan et al.l |200lt ); and 
due t o the much mor e prob lematic shear-ellipticity align- 
ment: IHirata fe Seliakl (|2004m ). Approaches that have been 
suggested so far to remove the contamination due to shear- 
ellipticity alignment suffer from several limitations and un- 
certaintie s: the statistical down-weighti ng of the effect pro- 
posed by Ijoachimi fc Schneider! |2008) relies on precision 
determinations of photometric galaxy reds hifts; other tech- 
niques, such as prop osedbve.g.lKingf<[2005l) . depend heavily 
on toy models (e.g. lCatelan et al j200lF which do not include 
any (as yet unknown) dependence of alignments on galaxy 
luminosity, mass, or redshift. Several studies using iV-body 
simulations however find evidence for a dependence of the 
stre ngth of halo-LSS alignment on mass and redshift (see 
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\20oS ) 

Unfortunately, it is currently still computationally un- 
feasible to use hydrodynamical simulations of galaxy forma- 
tion of the necessary resolution in a sufficiently large volume 
to investigate the correlation function of galaxy shapes and 
spins. To circumvent this barrier, some authors have used 



large-volume iV-body simulations together with a prescrip- 
tion to assign galactic disk orientations to the dark mat- 
ter haloes which typically either assumes perfect alignment 
with the DM halo spin or some random mis alignment angle 
between the two (e.g. iHevmans et al.| [2006). Also analytic 
models based on the halo model have been developed that 
include prescriptions for small s cale satellite-central align- 
ments (|Schneider fc Bridle! 12003 ) . Probability distributions 
for the angles between halo and galaxy spins have been de- 
termined by various authors (e.g. van den Bosch et alj|2002l; 
Baffin et al] 120051: Isharma fc Steinmet j 120051 : ICroft et all 
20081 : iBett et al.ll2009l ). ~ 

Accepting the current unfeasibility of large volume hy- 
drodynamic simulations, we propose a different approach. 
Combining local measurements of the alignment of galax- 
ies with the lar g e-scale structur e - as has been inves t igated 
by lHahn et all (|2007al lbh and lAragon-Calvo et alj (|2007h 
in pure iV-body simulations - circumvents the need for 
large volumes but can still provide constraints on the de- 
gree of galaxy align ment to be expected. The results of 
iNavarro et alj (12004! ) on the basis of four SPH simulations 
of the formation of single isolated galaxies indicate an align- 
ment of the spin of these galaxies with the intermediate prin- 
cipal axis of the tidal field at turnaround that persists to a 
slightly weaker degree also to z = 0. 

In this paper, we study the orientation of a sample of 
~ 100 galactic discs with respect to their larger scale envi- 
ronment in a hydrodynamical cosmological simulation of the 
formation of a large filament. We have chosen this partic- 
ular large-scale environment as it provides an environment 
where the gravitational shear is expected to be relatively 
large so that weak lensing shear should have good signal- 
to-noise. Any intrinsic alignment of galaxies on these scales 
will thus impact the measurements through the Gl-lensing 
term. Prior investigations in dark matter cosmological simu- 
lations have shown evidence for alignment of halo spin with 
filaments. However, considering the gas component rather 
than dark matter, the filament is a high pressure environ- 
ment so that hydrodynamic effects, such as ram-pressure, 
are expected to play a role in the evolution of the simulated 
galaxies. 

This paper is organised as follows: we describe the de- 
tails of our numerical simulation, the method to identify 
galaxies and the selection of galaxies for our analysis in Sec- 
tion [5] We first study the alignments between the gas, stel- 
lar, dark matter angular momenta and the total halo spin in 
Section [3] Then, we investigate in Section [J] the alignment 
of the various components with the cosmic large-scale struc- 
ture, quantitatively described by the tidal field. In SectionlAl 
we discuss the numerical convergence of our results, before 
a comparison of our results to observations is given in Sec- 
tion [5] Finally, we summarise and provide our conclusions 
in Section [5] 



2 NUMERICAL METHODS 

2.1 Specifics of the Cosmological Simulation 

We simulate the formation of a large cosmic filament and 
the galaxies em bedded in it using the AMR code Ramses 
(Teyssier 2002) to evolve dark matter, gas and stars. Gas dy- 
namics is computed using a second-order unsplit Godunov 
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Figure 1. The simulated filament at redshift z = 0. The images extend over 20/i~ 1 Mpc X 20/i _1 Mpc and have a projection depth of 
10 h~ 1 Mpc. (Left panel:) the dark matter particle distribution, (Right panel:) a composite RGB colour image of the gas distribution, the 
red image channel was assigned to the pressure, the green channel to the metallicity of the gas and the blue channel to the density. Since 
metal enriched gas is produced only in galaxies, each of the green blobs corresponds to a galaxy. The white circles in the right panel 
indicate the virial radii of those host haloes within the filament that have been excluded from refinement for most part of the simulation 
or are too close to a low-resolution region in the high resolution run (cf. Section 12.111 and are thus not considered in our analysis. All 
other galaxies are treated with full resolution and are not affected by low resolution outside of the zoom region. 



scheme, particles (dark matter and stars) are evolved using 
the particle-mesh method. The simulation includes s tandard 
recipes for star formation (R ascra fc Tevssierll2006l). super- 
nova feedback and chemical enrichment ( Dubois fc Tevssierl 
2008). We do not repeat details of these sub-grid models 
but kindly refer the reader to the mentioned papers. The 
filament was chosen by visual inspection of the environment 
of the most massive haloes in a dark-matter-only preflight 
simulation. We then chose a filament that is very prominent 
(and thus relatively massive) and not criss-crossed by too 
many smaller scale filaments. 

Initial conditions were generated using the Grafic-2 
tool from a power spectrum computed with Linger. We 
adopt cosmological pa rameters as obtained f rom the WMAP 
5th year data release jKomatsu et al.ll2009l ). i.e. we use the 
density parameters Q. m — 0.276, Slf, = 0.046, S1a = 0.724, a 
power spectrum normalisation of erg = 0.811 and a spectral 
index n 3 — 0.961 for the primordial spectrum, as well as 
a Hubble constant of Ho — 100 /ikms" 1 Mpc -1 with h = 
0.703. 

The filament is simulated using the "zoom-in" technique 
as a high-resolution region of 25/i _1 Mpc diameter at red- 
shift z — inside of a comoving 100ft _1 Mpc cosmological 
box. The effective initial resolution of our simulation is 1024 3 
particles and grid points inside this high-resolution region, 
dropping in steps of 8 to 128 3 in the remainder of the box. 
Additional refinements are confined to the high-resolution 
region only. We adopt a refinement strategy based on den- 
sity threshold so that the number of particles per cell (and 
the baryon mass per cell) remains roughly constant in time. 
Since the filament moves slightly in space, we introduce a 
"colour" field to mark the Lagrangian volume of the fila- 
ment. This field is passively advected during the formation 
of the structures and thus follows the evolution of the fila- 



ment perfectly. In order to save memory, we used as colour- 
ing variable the metallicity of the gas, with a very low (but 
non-zero) value of 10 -3 solar in the high resolution region 
and exactly zero outside. Additional refinements are only 
allowed where the metallicity exceeds one tenth of this ini- 
tial value. In addition, cooling and star formation are also 
switched off in non-coloured regions. The maximum resolu- 
tion of our simulation is physical 0.38 h~ 1 kpc at all times. In 
addition, our Lagrangian volume excludes the most massive 
galaxy groups (halo mass above 5 x 10 13 ft _1 MQ at z — 0) 
in order to speed up computation as we are only interested 
in smaller objects in this study. 

The two panels of Figure [1] give a visual impression 
of the filament at z = 0. We show two panels, the first 
illustrating the distribution of dark matter, the second is 
a composite image representing the density, pressure and 
metallicity distribution in the gas. Note that the main fila- 
ment is entirely wrapped by an accretion shock. It is thus 
clearly a high pressure environment. The less massive fila- 
ments remain cold except at node points where the more 
massive haloes reside. An inspection of the time-evolution 
of the simulation reveals that this accretion shock emerges 
first from the most massive haloes along the filament but 
then grows along the filament boundaries until the entire 
filament is enclosed by a shock. 

In the two panels of Figure [2] we show two examples of 
well-resolved disk-galaxies at z — in the simulation. The 
images represent a simulated stellar light composite image 
combining the i', r' and g' bands into an RGB image tak- 
ing into account dust reddening along the line of sight. The 
images also show the orientation of the stellar and gaseous 
angular momentum and indicate that these are perpendicu- 
lar to the disks and well aligned with one another. 
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centre of the gaseous disk is denned as the densest gas cell 
in the galaxy and the centre of the inner dark matter halo is 
defined as the most bound dark matter particle when eval- 
uating the potential for all dark matter particles that are 
contained in the galaxy patch. 

In order to determine halo masses and virial radii, we 
grow spheres around each galaxy centre until the enclosed 
matter overdensity falls below the overdensity A c of a viri- 
alised spherical perturb ation in a ACDM cosmolo gy as given 
by the fitting formula of lBrvan fc Normanl |l998). In a next 
step, we determine equivalence classes of galaxies that are 
contained within the virial radii of each other. The galaxy 
with the highest stellar mass in each equivalence class we call 
the "central galaxy", all other class members are "satellite 
galaxies" . We do not consider these satellite galaxies further 
in this study. 

Finally, unresolved galaxies are identified by tracking 
the dark matter particles that end up in galaxies back to the 
initial conditions and requiring that none of them started in 
a region that was not initially marked for cooling and re- 
finement. These unresolved galaxies are then removed from 
the galaxy catalogue. 

The central galaxies in our samples have stellar radii 
f"95 (which is the radius that contains 95 per cent of the 
stellar mass) that are at most 0.11 at z = 1, 0.095 at z = 
0.5 and 0.077 at z — of the virial radius of their parent 
halo. This makes us confident that the identified galaxies 
are not contaminated by other satellites in the halo and 
thus constitute clean samples of central galaxies. 




Figure 2. Simulated stellar light composite images of two well- 
resolved disk-galaxies using the i' , r' and g' band light and taking 
into account dust absorption. The white circles correspond to the 
radius containing 95 per cent of the stellar mass. The yellow line 
corresponds to the stellar angular momentum vector direction, 
the red line to the gas angular momentum. 



2.2 Identification of Galaxies 

In order to identify galaxies in the simulated data, we di- 
rectly use the AMR mesh structure and label contiguous vol- 
umes exceeding a density threshold. First, we pick the high- 
est refinement level which exceeds l/i _1 kpc resolution in or- 
der to smooth out galactic substructure. Next, we keep only 
those cells on this level, which exceed a threshold Sb = 50000 
in baryonic overdensity. This particular choice of St, has been 
found to be in excellent agreement with a visual identifica- 
tion of the galaxies. In a final step we label all contiguous 
regions that are left. This labelling is performed using a par- 
allel implemen tation of the Hoshen-Kopel man cluster find- 
ing algorithm l)Hoshen fc Kopelmanl [l97r3 l. Finally, all star 
particles and gas cells (also on higher refinement levels) are 
assigned to the mesh cluster Q^, i.e. the "galaxy", in which 
they reside. 

Galaxy stellar centres are then determined by itera- 
tively computing the stellar centre of mass within shrink- 
ing spherical apertures and moving the aperture to the new 
centre in each iteration. The procedure is stopped once less 
than 50 particles are contained within the aperture. The 



2.3 The Fiducial Galaxy Samples 

For our analysis, we keep only the best resolved galaxies with 
at least 25000 star particles, corresponding to a minimum 
stellar mass of ~ 7.5 x 10 9 fc _1 Mg. For these galaxies, rso, 
the radius containing 80% of the stellar mass is resolved by 
at least 8 resolution elements. Furthermore, this requirement 
implies that we can compare our results to a lower resolution 
run (see Appendix [A]) for which the minimum star particle 
number for galaxies of equal mass is then a factor of 8 lower, 
but still above 3000. 

We split the sample of well-resolved galaxies accord- 
ing to their host halo masses into 3 different halo groups. 
The high-mass sample is defined as all halos above the 
non-linear mass M*{z = 0) w 2 x 10 12 /i _1 M Q . Prior re- 
sults using pure JV-body simulations have indicated that 
M» might be a characteristic scale for the spin- alignment 
behaviour of dark matter haloes (see e.g. Hahn et al.ll2007bl ; 
lAragon-Calvo et al.l 120071 ; |Paz et aill200et l. The intermedi- 
ate sample is for halo masses smaller than this non-linear 
mass, but larger tha n the "cold accretion lim i t" of ~ 4 x 
lO 11 / ! " 1 ]^ fsee e.g. iBirnboim fc Dekell 120031 ; iKeres et al] 
120051 ; iDekel fc Birnboiml l200rj ). Haloes above this critical 
mass are believed to host a hot pressurized atmosphere that 
can exert ram pressure effects on the central galaxy. The 
third mass sample is for halo masses smaller than the criti- 
cal mass. Note that for the two higher redshifts z — 1 and 
z = 0.5, we merge together the high and intermediate mass 
samples to improve statistical significance. Detailed infor- 
mation about the galaxy samples is given in Table Q] The 
distribution of stellar masses in the samples is shown in Fig- 
ure [3] 
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LM MM HM 

z = 0: 

Afhalo / ft -1 Mo 1 - 4 x 10 11 5 x 10 11 - 2 x 10 12 > 2 x 10 12 

median M steUat / h -1 M s 1.3 x 10 10 6.8 x 10 10 4.2 x 10 11 

number of galaxies 70 39 12 

z = 0.5: 

A4aio / fe _1 M 1 - 4 x 10 11 - > 4 x 10 11 

median M stcllal . / fc^Mg 1.3 X 10 10 - 6.7 x 10 10 

number of galaxies 48 - 28 

z = l: 

Afhaic / ft~ 1 M Q 1 - 4 x 10 11 - > 4 x 10 11 

median M" stellar / ft _1 M 1.3 X 10 10 - 6.3 x 10 10 

number of galaxies 58 - 31 



Table 1. The galaxy samples used in our analysis. We split central galaxies according to their host halo mass in a low mass (LM) and 
high mass (HM) bin at redshifts 2 = 1 and z = 0.5 and an additional medium mass (MM) bin at z = 0. For each mass bin and each 
redshift, the median stellar mass as well as the number of galaxies in the bin is given. See Figure \3\ for histograms of stellar masses for 
each of the galaxy samples. 




Figure 3. Distributions of galaxy stellar masses in the two/three galaxy samples used in our analysis (see Table[T]for details) at redshifts 
2 = 1 (left), z = 0.5 (middle) and 2 = (right). 



2.4 Measured Galactic Properties 

In this study, we are interested in the orientation of disk 
galaxies. Since we exclude the most massive (large group- 
scale and above) haloes from refinement, indeed all galaxies 
considered are disky. We can thus simply quantify the disk 
orientation by the angular momenta of the various disk con- 
stituents: gas, stars and dark matter. 

We compute the angular momentum for each galaxy Qu 
by evaluating the sum 

Jfc = ^2 m,(r, - f fc ) x (v< - Vfc), (1) 

where ffc is the stellar centre of the fcth galaxy inner region 
(see previous section), and Vfc is the centre of mass velocity. 
The sum is evaluated separately for all star particles that 
are part of an identified galaxy C?fc to yield the stellar angu- 
lar momentum jf, as well as for all DM particles that are 
contained within a sphere of the radius of the most distant 
star particle in the galaxy J J . Finally, we also compute 
the gas angular momentum by summing over all leaf 
cells of the AMR tree that are contained in the galaxy (and 
are thus effectively bounded by an iso-density surface - see 



Section \2. 2 1 for details). For gas, ffc refers to the densest cell 
in the galaxy patch and Vfc to the centre of mass velocity 
computed over all leaf cells in Qk, while for dark matter ffc 
refers to the most bound dark matter particle and vfc to the 
centre of mass velocity of the dark matter particles in the 
sphere containing (see also Section [2. 2 \ . In order to de- 
termine the total halo angular momentum, we evaluate eq. 
JT} summing over all components and up to the virial radius 
to obtain Jtot- 



2.5 Quantifying Large-Scale Structure with the 
Tidal Field 

In ana logy to the approach followed by Irlahn et al.l 
i|2007al lbl). we use the eigenvectors of the tidal field ten- 
sor smoothed on some scale R 3 to probe alignments with 
the large-scale structure. This app roach has t w o mai n 
advantages: (1) as demonstrated in lHahn et alj (|2007al ). 
the tidal field provides a natural dynamical quantifica- 
tion of the geometry of the large-scale structure; and 
(2) any local alignments with the tidal field directly 
translate into shear-ellipticity alignments in weak lensing 
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Figure 4. Directions of the tidal field eigenvectors at z = 0.5 for a smoothing scale of R s = 2?t — 1 Mpc physical (3/i _1 Mpc comoving). 
The V3 eigenvector indicates the overall direction of the filaments, while V2 and vi span the perpendicular plane. Linear tidal torque 
theory favours alignment of galactic spin with V2 . 



ijHirata fc Seliakl |2004 >. In addition, when supplemented 
wit h auto-correlation functions of the non-linear tidal field 
(see iLee et al l I2009al 1 , local alignments with the tidal field 
can be translated into correlation functions. We will inves- 
tigate this approach in future work. 

We determine the tidal field by first constructing the 
full matter overdensity field <5(x) of the entire simulation 
volume on a 512 3 grid. This grid is subsequently smoothed 
with a Gaussian kernel of varying filtering radius R a . Finally, 
we solve Poisson's equation using the fast Fourier transform 
(FFT) to obtain the potential 0(x), evaluate the Hessian 
matrix Tij = didj(f> by finite differencing the potential and 
interpolating via cloud-in-cell to the galaxy positions. Di- 
agonalising Tij at each galaxy's position, we find the three 
eigenvalues \i(R a ) < \2{R a ) ^ \\{R S ) and the correspond- 
ing unit eigenvectors V3(i? s ), V2(-R s ), vi(f? s ) as a function 
of the applied filtering scale R„. Note that we define Tij to 
be simply the Hessian of the potential, i.e. the signs of the 
eigenvalues might be reversed with respect to some other 
studies. 

The eigenstructure of the Hessian of the gravitational 
potential is intimately connected to the large-scale dynam- 
ics of the cosmic web. Filamentary regions tend to have one 
negative and two positive eigenvalues when smoothed on 
large-enough scales. In particular, we find e.g. average values 
(A3} = —0.2, (A2) = 0.2 and (Ai) = 0.6 on scales of comov- 
ing 3 h~ 1 Mpc where the averages are taken over the sample 
of our galaxies at z = 0.5. Units are such that Xi = 5, 
where S is the matter overdensity. The signature of these 
eigenvalues indicates that the galaxies are, on average, in- 
deed embedded in a filamentary environment. The signs of 
the eigenvalues reflect the large-scale expansion along the 
filament - corresponding to the direction indicated by the 
eigenvector of the negative eigenvalue A3 - and gravitational 
contraction in the perpendicular plane - corresponding to 
the plane spanned by the eigenvectors V2 and Vi. The con- 
nection between the eigenstructure of Tij and dynamics fol- 
lows readily from the Zel'dovich theory of pancake collapse 
(|Zel'Dovichlll970l ). In Figure U we show the orientations of 
the eigenvectors with respect to the cosmic web structure. 

To summarize, V3 indicates the overall direction of a 
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Table 2. (Upper half) Median alignment angles between the gas 
(G), stellar (S) and dark matter (DM) component of the galactic 
disk at redhifts 2 = 0, 0.5 and 1. (Lower half) Median alignment 
angles between the three components of the galactic disk and 
the total angular momentum of the host halo. Errors are non- 
symmetric uncertainties on the median (see text for more details). 
Values are given for the entire sample of well-resolved central 
galaxies since no mass dependence was found. 



filament, while V2 and vi span t he perpendicular pl ane. A 
more detailed account is given in lHahn et alj |2007al ). 



3 THE DISTRIBUTION OF ANGULAR 
MOMENTUM IN GALAXIES 

3.1 Angular momentum correlations within the 
galaxy 

In Figure O we show the distributions of angles 9 between 
the angular momentum vectors of the stellar (S), gas (G) and 
dark matter (DM) component within well-resolved central 
galaxies at z = (the nomenclature is such that, e.g., #s,g 
is the angle between the stellar and gas angular momentum 
vectors). The corresponding median angles are given for 2 = 
— 1 in Table [2] (upper half) together with the error on the 
median. We compute the upper error a u and lower error ai 
as 
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Figure 5. Distribution of angles between the angular momentum 
vectors of the stellar, gas and dark matter component within the 
blob that identifies the galaxy. The angle distributions are shown 
for the entire sample of well-resolved central galaxies. 



iBett et all (|2009h who report a smaller median angle of 34 
degrees between galaxy and outer dark matter halo at z — 0. 
It remains a matter of future studies to better quantify the 
influence of resolution and feedback recipes on these results. 
A further possible origin for the disc repancy between our 
results and those of IBett et al.l (|2009T ) might lie in the dif- 
ference of large scale environment. While we simulate the 
relatively high-pressure environme nt of a mas s ive fi lament, 
the galaxy sample considered by IBett et ail (120091 ) has a 
halo mass upper limit of 1.6 x 1O 12 M0, suggesting that they 
reside in a lower density environment with a more quiescent 
formation history. 

The outer and inner dark matter spins are significantly 
stronger aligned than the baryonic components. We find a 
median angle of ~ 45 degrees at z — which decreases 
to ~ 34 degrees at z = 1. Again, we find no correlation 
between the alignment angles and either stellar or host halo 
mass with a significance above la. 



where subscripted numbers indicate percentiles and N is the 
number of galaxies considered. 

Within the galactic disks (which is defined by the mesh 
cluster identified as a single galaxy - see Section [2.21 - and 
thus has no imposed symmetry for the stars and the gas), the 
gas angular momentum deviates by a median of 8 degrees 
from that of the stars. The difference with respect to the 
dark matter angular momentum, which was determined in 
a sphere around the disk, is slightly higher at around 18 
degrees at z = for both stars and gas and increases to 
around 36 degrees at z = 1. 

We find no correlation between these alignment angles 
and either stellar or host halo mass with a significance above 
la. This is a comforting result as stellar mass is directly 
proportional to the number of star particles and the lack of 
a correlation demonstrates the absence of an obvious res- 
olution problem. We discuss numerical convergence of our 
results in more detail in Appendix IA1 



3.2 Angular momentum correlations between 
central galaxies and their host halo 

In Figure |S] we show the distribution of angles between the 
angular momentum vector of the host dark matter halo (de- 
termined out to the virial radius from both baryonic and 
dark matter) and the inner components (determined at the 
disk) for the three redshifts at which we perform our anal- 
ysis. Galaxies from all samples have been combined. The 
corresponding median angles and the width of the distribu- 
tions are given in Table [2] 

We find a fairly broad distribution of angles with the 
host angular momentum deviating by ~ 50 degrees at z — 
from either the stellar or gas spin of the central galaxy. Our 
results indicate that this discrepancy decreases slightly to 
about 46 degrees at z = 1. The value we find for the stellar 
co mponent is in excel lent agreement with the results found 
by ICroft et al.l (|2008h who report a median angle for the 
stellar component of 44 degrees and 70 degrees for the gas. 
The quoted larger angle for the gaseous component is pos- 
sibly due t o sate llites that are included in the sample of 
ICroft et al.l (|2008l ) and for which ram pressure modifies the 
gas disk orientations very quickly. It appears, however, that 
our results are in slight tension with the recent findings of 



4 THE LARGE-SCALE ORIENTATIONS OF 
GALAXIES 

4.1 Predictions of Analytic Models 

The linear tidal torque theory predicts that the angular mo- 
mentum of a proto-galaxy is determined by its geometric 
shape and the tidal field exerted by the larger scale environ - 
ment (e.g. |Peebles|[l969l ; IWhitelll984l ; iPorciani et ai1l2002r i. 
In particular, angular momentum is generated from the mis- 
alignment of the Lagrangian moment of inertia tensor Iij 
and the Lagrangian tidal tensor Tj, so that at first order 
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where is the Levi-Civita symbol. In the eigenspace of 
the tidal tensor, the components of J become 



(4) 



Since, by definition, A3 — Ai is the largest coefficient, lin- 
ear tidal torque theory thus predicts that, at first order, J 
is preferentially aligned with V2 in a statistical sample of 
haloes/galaxies where contributions from the moment of in- 
ertia tensor average out and Tij and Jy are assumed to be 
uncorrelated. Moreover, generation of angular momentum is 
most efficient at turnaround when the Iij is maximal. 

However, the assumption that Tij and Iij are uncorre- 
lated is not well justified, as the two tensors are indeed found 
to be correlated (see e.g. IPorciani et al.l |2002| ; iLee et all 
2009b). Thus it is particularly interesting to see whether 
the alignment predicted from Lagrangian perturbation the- 
ory is indeed present in hydrodynamic cosmological simula- 
tions, and particularly so at late times well after turnaround 
of t he proto-galaxy. 

ILee fc Erdogdul (|2007h have derived a probability dis- 
tribution for the angle cos #2 = |J • V2|/||J|| - along which 
the possibly strongest alignment could be observed - of the 
form 



p(cosr? 2 ) = (l + c) A /l-- 



I — 1 • ( 1 — — cos 2 82 



-3/2 



(5) 
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Figure 6. Distribution of angles between the total (dark matter plus baryonic) angular momentum vectors of the dark matter host 
halo and the stellar, gas and dark matter component of the central galaxy. The angle distributions are shown for the entire sample of 
well-resolved central galaxies at z = (left), z = 0.5 (centre) and 2 = 1 (right panel). 



where c is a normalisation parameter. 

4.2 The Distribution of Spin and Tidal Field 
Eigenvectors 

4-2.1 Anisotropy of the Tidal Field 

In Figure [7] we show the spatial distribution at z = 0.5 
of tidal field eigenvectors V1...3 for the comoving smoothing 
scales R s — 1,2 and 5/i _1 Mpc in Mollweide projection. The 
tidal field is evaluated at the positions of the galaxies that 
are included in any of the samples (LM or HM). As this is an 
equal-area projection, local concentrations directly translate 
into an increased anisotropy for that direction. Since the 
tidal field is a tensor field, each vector has spin-2 symmetry 
and hence is represented by two points, the second being a 
reflection of the first at the origin. While the tidal field is 
already highly anisotropic on scales of 1 /i -1 Mpc, on scales of 
5/i _1 Mpc virtually only the main filament remains so that 
the distribution of eigenvectors becomes almost dipolar. 

4-2.2 Anisotropy of the Galaxy Spin Distribution 

Figure [8] shows the distribution of stellar angular momentum 
vectors at redshifts z = 1, 0.5 and in Mollweide projection. 
At all redshifts the distributions are clearly anisotropic. The 
origin and detailed form of this anisotropic distribution is 
determined by both physical and numerical effects. We will 
discuss their interaction in what follows. 

At z = 1, the distribution is anisotropic with no clear 
clustering of spin vectors around the Cartesian grid vectors. 
This provides strong evidence that the high-redshift spin 
distribution is anisotropic due to physical processes. Com- 
paring with Figure 6, the tidal field on scales < 5/i -1 Mpc 
is anisotropic at a similar level. We will discuss their cross- 
correlation in the next section. 

At z = 0.5 there is a slight tendency for the spins 
to cluster around the grid vectors that is more obvious 
at 2 = for the low-mass galaxies. Since a multi-grid 
Poisson solver (such as the one used in Ramses) intro- 
duces preferred directions along the Cartesian grid ba- 
sis vectors, a non-physical anisotropy of the gravitational 
inter-particle force is expected close to the resolution limit 



(cf. Irlocknev fc Eastwood! Il98ll ) arising from an unavoid- 
able lack of symmetry in the discrete Laplacian operator. 
These errors are sy stematic in any asymmetric system (cf. 
iMav fc Jameslll984 ). This is clearly a systematic effect in our 
numerical results, and we need to estimate its magnitude. 

Remarkably, the grid-aligned galaxies at low redshift 
are not distributed isotropically among all Cartesian grid 
vectors. In larger scale simulations of entire cosmological 
volumes (i.e. not restricted to a fila mentary region - as e.g. 
in the "Mare Nostrum" simulation. lOcvirk et al.ll200ct ). we 
do not see this effect. In these larger volume simulations at 
lower spatial resolution, grid-alignment is also observed but 
with galaxies aligned with all three grid vectors. 

Given the highly anisotropic distribution of tidal field 
eigenvectors (see previous subsection), it is however clear 
how such an anisotropy can arise. Assuming that the galax- 
ies' orientations are determined by some anisotropic physi- 
cal effect (such as tides), the disks' initial angular momen- 
tum will be anisotropic (as it is at z = 1). Given the ad- 
ditional numerical anisotropy, the disks will subsequently 
relax slowly to the closest potential minimum, i.e. the clos- 
est Cartesian grid vector, in the absence of any perturba- 
tion that is larger than the numerical error (and thus pre- 
dominantly in low density regions). This will happen on a 
timescale that is connected to the anisotropy of the error 
in the multi-grid Poisson solver. Support for this interpre- 
tation of the numerical results comes from a look at the 
galaxies' spin distribution at high redshift where no obvious 
grid-alignment is observed. Looking at the left panel of Fig- 
ure [5] we see a non-isotropic distribution with no spins in 
the vicinity of those grid axes where by z = also no grid- 
aligned galaxies are seen (e.g. the one in the very centre of 
the Mollweide maps). 

We will continue to discuss the impact of such a spu- 
rious alignment on our results in Section [X] To conclude 
this section, we can however clearly say that the distribu- 
tion of galaxy angular momenta is anisotropic due to galaxy 
formation physics. The anisotropy is then likely amplified 
by numerical effects at low redshifts, but numerical effects 
alone will never generate an anisotropic distribution. 
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Figure 7. Mollweide sphere projections of the three eigenvalues of the tidal field at the positions of the galaxies for three smoothing 
scales, lh —1 Mpc (top), 2h — 1 Mpc (middle) and 5/i — x Mpc (bottom). All data is for redshift z = 0.5. Since the tidal field has a spin-2 
symmetry, each vector is represented by two points. The Mollweide projection is an area preserving projection. 




Figure 8. The orientations of galaxy stellar angular momenta at redshifts 2 = 1 (left), z = 0.5 (middle) and 2 = (right). The angular 
momentum vectors are shown in the Mollweide projection. The circled crosses indicate the Cartesian basis vectors of the AMR grid. 



4.3 Correlations with the tidal field 

As outlined in Section 14.11 linear tidal torque theory indi- 
cates that galaxies might show a preferential alignment with 
the intermediate principal axis of the tidal field. In this Sec- 
tion, we investigate whether such an alignment persists also 
at times after the turnaround of the proto-galaxy in our 
simulation. Unfortunately, the scale on which the tidal field 
should be computed is not known a priori so that we present 
the degree of alignment for a range of smoothing scales. 
Since these scales are likely to be different for galaxies of 
different masses, and galaxies of different halo masses may 
also experiences differences in the details of their formation, 
all results are split into different halo masses as described in 
Section secigalsamples. In addition, we distinguish low-mass 
galaxies in high-density and in low-density environments. 

In Figure [5] we show the scale dependence of the cosine 
of the angles | cos 9i | between the angular momentum vec- 
tors of the galaxies and the tidal field eigenvectors for the 
inner dark matter, gas and stellar components. The panels 
represent the alignments with the 3 eigenvectors (vertical) 
for the 2 mass bins (horizontal). We additionally split the 
low mass bin (LM) into an overdense bin 8 > <5so and an un- 
derdense bin <5 < S^o, where here 5 is the matter overdensity 
field smoothed on a 2/i _1 Mpc physical scale and Sso is the 



median overdensity for the galaxies in the respective mass 
bin. The values of <5so for our sample of galaxies turn out to 
be: 0.27 at z = 1, 0.54 atz = 0.5 and 1.16 at z = 0. 

Using a two-tailed Kolmogorov-Smirnov test, we as- 
sess the inconsistency of the angle distributions on a scale 
of physical 2/i _1 Mpc with a fiat (i.e. random distribution) 
Prand(| cos 8i | ) = 1. The significance levels are given in units 
of o in Table El 

At 2 = 1, we find the strongest alignment for the 
low-mass low-density sample, where on scales of physical 
2/i~ 1 Mpc, the alignment signal with the intermediate prin- 
cipal axes of the tidal field peaks at a median angle of 
j cos 02 1 ~ 0.7 for the stellar angular momentum component. 
The deviation from a random angle distribution is significant 
at 2.75<r. The alignment of the gas disk is slightly weaker, 
while that of the inner dark matter is slightly stronger. There 
is a weaker counter-alignment with the third principal axis 
on these scales, while the alignment with the first axis is 
consistent with random. 

For the low-mass high-density sample, we find a sig- 
nificant alignment only with the first principal axis for the 
stellar and gas component, while the dark matter angular 
momentum is consistent with random at almost all scales. 

Finally, for the high-mass sample, we find a significant 
alignment of the stellar and gas component with the third 
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Figure 9. Median cosine at z = 1 of the angle between the spin of the dark matter (red), gaseous (green) and stellar (blue) component 
of central galaxies and the three eigenvectors of the tidal field tensor as a function of the scale R s (given in h^ 1 Mpc comoving) on which 
the tidal field was smoothed. The data is shown in for the halo mass bins given in Table [l] with the low mass (LM) sample split into 
two samples above and below the median environmental matter overdensity determined on physical scales of 2/i _1 Mpc (i.e. 4/i _1 Mpc 
comoving). The shaded areas represent the la uncertainty on the median obtained with 200 bootstrap resamplings. 
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Table 3. Significance levels for the rejection of the null-hypothesis that the angle distributions between the spin of the gas disk (G), 
stellar disk (S) or inner dark matter (DM) with the tidal field eigenvectors are consistent with a uniform distribution p(| cos6>;|) = 1 
obtained with the Kolmogorov-Smirnov test. The tidal field was smoothed on a scale of physical 2 ?t~ 1 Mpc for this test. Boldface numbers 
highlight those correlations that have a significance above 2a. 
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principal axis on comoving scales of 3 ft - Mpc which is how- 
ever not seen for the inner dark matter component at z — 1. 
At lower redshifts, this alignment decreases slightly but is 
also seen for the dark matter component. 

In Figure [TOl we show the scale dependent median align- 
ment for the z = 0.5 samples of galaxies. As for the 2 = 1 
sample, we again find a preferential alignment of the low- 
mass low-density sample with the intermediate axis of the 
tidal field peaking on scales of physical 2ft _1 Mpc and a 
counter-alignment with the first principal axis. For the low- 
mass high-density environment, our findings are consistent 
with random. There is again a weak tendency for the high 
mass sample to be aligned with the third principal axis. 

In Figure rm we show the scale dependent median align- 
ment for the z = samples of galaxies. For the low-mass 
samples, all alignment signals are consistent with random. 
For the medium mass sample, we find a very weak ten- 
dency for alignment with the intermediate principal axis 
of the tidal field. The high-mass sample shows a signifi- 
cant alignment with the third principal axis on small scales 
< 1 /i~ 1 Mpc that is still present but weakens at larger scales. 
Note however, that the high-mass sample consists only of 12 
galaxies. 

4.4 Alignment Angle Probability Distributions 

Figures 1121 1131 and Q3] show the probability distributions 
p(|cos#i|) for the low mass samples (LM) at redshifts z = 
1, 0.5 and 0, respectively. The galaxy samples have again 
been split into a high and low environment density sample 
above and below the median density, as described in Section 
14.31 and the angle distributions are given for the same scale 
on which the density field has been obtained - 2/i _1 Mpc 
physical at all redshifts. 

For the higher redshifts z — 1 and 0.5, the low den- 
sity sample is consistent with the model fit for spin align- 
ment with the interm ediate axis of the tidal field given by 
iLee fc Erdogdul |2007l ) albeit for a parameter which is sig- 
nificantly larger (roughly a factor of 9) than the best fit 
value to the 2MASS galaxy samples obtained in their study. 
While this discrepancy might seem large, it is not clear to 
what extent our samples are comparable at all with the ob- 
servational data from the 2MASS survey. Furthermore, the 
weaker alignment seen in observational data might be due 
to the inherent difficulty to infer the orientation of a galaxy 
from the projected image, as well as uncertainties in deter- 
mining the tidal field. 

More curiously however. ILee fe Erdogd u (2007) find a 
stronger alignment signal with the intermediate axis for 
galaxies in high density regions. At all redshifts, for the high- 
density sample of low-mass galaxies, our results are consis- 
tent with random orientations. While tidal torque theory 
predicts that alignment with the intermed iate axis should 
indee d be stronger in high density regions l|Lee fc Erdogdul 
l2007h . our results do not support that this prediction is valid 
in highly non-linear environments. 

4.5 The Impact of Environment on Alignments 

In Figure [T5l we show the dependence on environmental den- 
sity of the alignment angles of the galaxies' stellar angular 



momentum with the three principal axes of the tidal field for 
the low- mass sample of galaxies at redshift z = 1. The over- 
density field S is evaluated on the scale of R s = 2/i _1 Mpc 
physical. 

We find that the alignment of stellar angular momen- 
tum with the intermediate axis decreases strongly when the 
environmental density increases. This is perfectly consistent 
with our the results from Figure [9] that no alignment with 
the intermediate axis is present for haloes in regions above 
the median environmental overdensity. On the other hand, 
the alignments with the first and the third axis increase with 
the overdensity of the environment. 

All three relations between 5 and the three | cos Qi | are 
significant. In particular, we find Spearman rank correlation 
coefficients of 0.32 with 2.47<r significance for 9$, —0.38 with 
2.96a significance for 02 and 0.24 with 1.81<r significance for 
61. We find quantitatively identical but weaker and slightly 
less significant correlations also at z — 0.5. 

These results strongly suggest that low-mass galaxies 
are initially aligned with the intermediate principal axis, as 
predicted by tidal torque theory, and that this alignment is 
subsequently destroyed, possibly by environmental effects. 
Since the alignment with the third principal axis increases 
during this process, it is plausible that anisotropic galaxy 
mergers in the filaments or the galaxy's motion along the 
filament and associated ram-pressure effects or a tidal in- 
fluence on the ga s flow accreting onto the galaxy (see also 
lHahn et alj|200a ) reorient the angular momentum. A first 
analysis of correlations of the galaxies' mass accretion rates 
(as a proxy for mergers) and ram-pressure with alignment 
strength did not lead to conclusive results. 

Since tidal torque theory predicts that angular mo- 
mentum is generated most efficiently at turnaround, it is 
also plausible that the galaxy's motion simply leads to a 
decorrelation of the tidal field it resides in at later times 
from that it experienced at turnaround. Galaxies in higher 
density regions are likely to have moved further. Specifi- 
cally, the low-mass galaxies at z = 1 have moved from 
their Lagrangian position (i.e. their centre of mass posi- 
tion in the initial conditions of the simulation) on average 
3.3±1.2 ft _1 Mpc, comoving, when they reside in regions with 
8 < S50 , while the overdense sample shows a mean displace- 
ment of 3.9 ± 1.2 ft Mpc. We find correlation coefficients 
between displacement and 62 of —0.25. Given that there is 
no ambiguity in smoothing scale when using the displace- 
ments, it is to be expected that this correlation coefficient 
should be much larger if decorrelation were the only reason 
for the deviation from tidal-torque prediction for the high- 
density sample. 

Given the inconclusive evidence at this stage, we post- 
pone a thorough investigation of the origin of the observed 
density dependence to future work. 



5 COMPARISON WITH OBSERVATIONS 

Most observational results for the alignments of central 
galaxies exist for elliptical galaxies. Early studies found an 
alignment of brig htest cluster ga laxies with the shape of the 
cluster itself (e.g. lBinggeli|[l982T ). More recently, the align- 
ment of luminous red galaxies (LRGs) has been detected 
at high significance in the Sloan Digital Sky Survey(SDSS) 
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Figure 10. Same as Figure[9]but at redshift z = 0.5. 



<|Mandelbaum et al.ll2006l ; iHirata et al.l 120071 ). These LRGs 
show radial alignment with their surrounding density peaks. 
These observational detections of elliptical galaxy align- 
ment are consistent with theoretical studies that demon- 
strate that they form through anisotropic mergers alo ng the 
filaments at whose n odes they will eventually reside |Westl 
1 19891 ; |Pubinskilll998T ). Our results indicate that the spin of 
the most massive disk galaxies in our samples show a pref- 
erence for alignment with the third principal axis of the 
tidal field which is parallel to the overall direction of the 
filaments (cf. Figures [9l 1 101 ITT|) . This thus corresponds to 
a radial anti-alignment of the disk shape with the filament. 
Since we however excluded the scales of galaxy groups from 
our simulations, we can only speculate on their alignment 
properties. It is nonetheless perfectly consistent with obser- 
vational results that the anti-aligned disk galaxies will merge 
along the massive filament with the merger remnant then 
being aligned with the filament since the orbital angular 
momentum of a major merger will outweigh the contribu- 
tion from the spin a ngular momentu m of the contributing 
disk galaxies (see also |Pubinski|[l99a ) . Our result could also 
ex plain why the alignment s ignal of massive ga l axies seen 
by iMandelbaum et all |200fJ ) and IHirata et all (|2007r ) de- 
cays so rapidly when less luminous galaxies are considered 
and thus the fraction of massive disk galaxies increases. 

For lower mass disk galaxies, the observational find- 



ings a re h owever less c o nclus ive. Both Mand elbaum et al.l 
( 2006) and IHirata et all (|2007n detect no al ignment for L < 
L* ga laxies in SDSS. On the other hand, iLee fc Erdogdul 
(2007) claim a detection of an alignment of 2MASS galax- 
ies with the intermediate principal axis of the recon- 
structed tidal field. Such an alignment with the interme- 
diate axis is consistent with an alignment of galaxy rota- 
tion axes with the shells around cosmic voids. The detec- 
tion of such an alig nment of disk galax ies around voids 
has been claimed by iTruiillo et all i|2006h in the 2dFGRS 
and SDSS galaxy surveys. Their analysis was recently re- 
peated by ISlosar fc White] (|2009l) on the SPSS D R6 who 
how ever cannot repr o duce the Truiillo et al.l (|2006l ) result. 
The [1 osar fc White! (2009) analysis extends out to a red- 
shift of < 0.18. The strong density dependence of disk galaxy 
alignment that we found in our analysis is likely to impact 
these analyses of volume limited galaxy samples. We thus be- 
lieve that our results are consistent with their non-detection. 
It is also plausible that the alignment is weaker at low red- 
shifts. 



6 SUMMARY AND CONCLUSIONS 

We have used a cosmological hydrodynamic simulation of 
galaxy formation in a large-scale filament to investigate (1) 
the orientations of the angular momentum of the stellar and 
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Figure 11. Same as Figure|9]but at redshift 



gaseous disk of central galaxies, as well as the inner part of 
their dark matter halo with respect to the total angular mo- 
mentum of the host halo; and (2) the orientations of their 
angular momentum with the surrounding large-scale struc- 
ture by computing scale-dependent alignments with the tidal 
field eigenvectors. 

The simulation provides us with a sample of ~ 100 
disk galaxies at 380 h~ pc resolution down to z = 0, span- 
ning halo masses between 10 11 and 10 13 li _1 Mg and stellar 
masses between 7.5 x 10 9 and 8 x 10 11 ft _1 M . We focus this 
first analysis on the three simulation snapshots at redshifts 
0, 0.5 and 1. We split our sample of galaxies into low-mass 
galaxies, which are those with a halo mass below the cold 
accretion limit of ~ 4 x 10 11 /i _1 Mq, and high-mass galaxies 
above this limit. At z — 0, we consider an additional medium 
mass bin bridging the gap between the cold accretion limit 
and the non-linear mass scale of ~ 2 x 10 12 /i _1 Mq. Our 
main results are summarised as follows: 



• There is an almost perfect alignment at a median of 
~ 18 degrees of the stellar, gaseous and inner dark mat- 
ter angular momenta at low redshifts. At z = 1, there is a 
slightly weaker alignment at ~ 36 degrees of the stellar and 
gaseous spins with the dark matter spin, likely due to the 
higher fraction of unrelaxed galaxies at that epoch. We do 



not find any dependence of this alignment signal either on 
environmental density, or halo mass or stellar mass. 

• The distribution of angles between the spin of the cen- 
tral galaxy and the entire host halo is significantly broader, 
the corresponding median angles larger. We find a median 
angle of ~ 50 degrees between both the stellar and the gas 
disk and the total halo angular momentum at z = 0. This 
median angle decreases slightly at higher redshifts to a me- 
dian of ~ 46 degrees at z = 1. The spin of the inner dark 
matter halo is slightly stronger aligned with the total halo. 
We find a median of ~ 45 degrees at z = and ~ 34 degrees 
at z = 1. Again, there is no dependence of this alignment 
signal on either stellar or halo mass, or on environmental 
density. 

• Low-mass galaxies in low density regions are aligned 
with the intermediate axis of the large-scale tidal field ten- 
sor, peaking on scales of physical 2/i _1 Mpc (i.e. 4/t -1 Mpc 
comoving) at z = 1 with very high significance. Such an 
alignment is consistent with the predictions from linear tidal 
torque theory. 

• In density regions above the median overdensity, on 
scales of physical 2/i _1 Mpc, we find however no evidence 
for alignment of the low-mass galaxies with the intermedi- 
ate principal axis of the tidal field - at any redshift since 
z = 1. Instead, we find alignment with the first principal 
axis at z = 1, which disappears at later epochs. 
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Figure 12. The distribution of | cos6\| between the galaxy spin vector and the three eigenvectors V3 (left), V2 (middle) and vi (right) 
at z = 1. Shown is the binned median, errorbars correspond to the Poisson error on the median. The distributions are shown for the 
stellar (blue circles), gaseous (green squares) and dark mat ter (red diamonds) components of the galaxies. The dashed and dotted line 
indicate the mod el predictions from | Lee fc Erdogdul i2007T ) .The dashed line corresponds to the observationally determined fit of the 
model as given in I Lee fc Erdogdul 1 120071 ). while the dotted line corresponds to the model with c = 0.7. In each plot, the top panel shows 
the distribution for the low-mass low-density (LM, <5 < (S50) sample, and the lower the distribution for the low-mass high-density (LM, 
<5 > 55o) sample. Density and alignments have been determined on scales of 2h _1 Mpc physical. 
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Figure 13. Same as Figure [T2l but for the low mass (LM) samples at 2 = 0.5, split in low (top) and high (bottom) environmental density. 



• We find a strong and significant correlation between 
environmental density and the degree of alignment with the 
intermediate axis at z = 1. Low mass galaxies gradually 
misalign with the intermediate principal axis of the tidal 
field with increasing environmental density. A possible ori- 
gin of this effect could be either torques exerted by ram- 
pressure when the galaxies enter a high-pressure environ- 
ment, a change in accretion/merger rates in these regions, a 
gradual decorrelation from the tidal field at formation due 
to the galaxy's motion in dense regions, or a combination of 
these effects. We will quantitatively investigate this issue in 
a future study. 



• The alignment with the intermediate principal axis of 
low-mass galaxies in low-density regions weakens with de- 
creasing redshift. The interpretation of this result is not 
straightforward, as low mass galaxies show some spurious 
alignment with the AMR mesh at low redshifts. It is im- 
portant to stress however that we have evidence, from the 
analysis of the z = 1 and z = 0.5 snapshots, that non- linear 
effects do weaken the correlation, raising the question as to 
whether the absence of alignment at z = is a spurious 
effect due to numerical effects in our simulations, or rather 
an important physical effect. Simulations are planned which 
are better tailored to establish whether the alignment with 
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Figure 14. Same as Figure 1121 but for the low mass (LM) samples at z = 0, split in low (top) and high (bottom) environmental density. 
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Figure 15. Cosine of the alignment angle of the stellar angular momentum with the three principal axes of the tidal field for all galaxies 
in the LM sample at z = 1 as a function of environmental density determined on scales of physical 2/i -1 Mpc. Light gray points indicate 
the data, dark gray points show binned median values with errorbars representing the uncertainty on the median. The lines represent 
linear least square fits to the data. 



the intermediate axis persists in low density regions down 
to z = 0. Simulations of isolated galaxies hint at the possi- 
bility that an alignm ent signal could persist down to z — 
l|Navarro et al.ll2004r ). 

• Within the caveats of a limited statistical significance, 
due to small number statistics, an important trend is un- 
veiled, for the most massive disk galaxies at all redshifts to 
be aligned with the third principal axis of the tidal field 
(pointing along the filaments). 



Summarising, low-mass galaxies with halo mass below 
4 x 10 11 h~ 1 Mpc are likely to show strong alignment with 
the intermediate principal axis of the tidal field at redshifts 
relevant for weak lensing studies. This alignment is signif- 
icantly weakened in regions of high environmental density, 
likely due to non-linear effects. We are studying the implica- 
tions of such alignments for weak lensing experiments such 
as the ESA mission Euclid, under study for launch before 
the end of the decade to start the quest to constrain the 
nature of Dark Energy. 
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APPENDIX A: ASSESSMENT OF NUMERICAL 
EFFECTS ON GALAXY ORIENTATIONS 

Al The Lower Resolution Simulation 

We also performed a simulation of the filament described 
in Section 12.11 at identical mass resolution but with a lower 
maximum refinement level so that the maximum spatial res- 
olution is 0.76/i~ 1 kpc at all times, a factor of 2 lower than 
the simulation used in our analysis. We analysed this simu- 
lation in an identical way as the higher resolution one and 
quantify the differences in what follows in order to deter- 
mine the numerical reliability of our results. The lower spa- 
tial resolution should especially demonstrate whether the 
galaxy orientations are affected by force anisotropy at the 
resolution level. 

We decide to keep the mass samples identical to those 
in the previous sections so that galaxies on average contain 
8 times less star particles (and thus more than 3000) but 
an equal amount of dark matter particles. We compare our 
results at z — as spurious grid alignment effects are ex- 
pected to be strongest in the dark energy dominated period 
where the merger rates decline substantially and the galaxies 
thus could potentially relax to the grid. The low-resolution 
simulation contains slightly more galaxies since additional 
refinement, cooling and star formation were not restricted 
to not contain the massive groups. These group galaxies are 
however not included in our comparison as we restrict our- 
selves to the central galaxies of haloes that are resolved in 
both simulations. 

A2 Angular Momentum Correlations in the 
Galaxy and with the Host Halo 

In Table IAU we give the median angles between the stel- 
lar, gas and dark matter component within the galaxies and 
with the total halo angular momentum at z — for the 
low-resolution run. Comparing these angles with those of 
the high-resolution run in Table [2l we find that the me- 
dian alignment angles obtained at lower spatial resolution 
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Table Al. Same as Table[2]but for the simulation with two-times 
lower spatial resolution. 




Figure Al. Same as Figure [8] but for the low-resolution run at 

2 = 0. 

are marginally larger, the corresponding vectors thus less 
aligned. Given the width of the distributions, the results 
can however be considered consistent at high significance. 

A3 Isotropy of the Galaxy Spins 

In Figure IAU we show the orientations of stellar angular 
momentum vectors for the low-resolution simulation. This 
should be compared with the corresponding Figure [8] for the 
high-resolution run. Visually, we see a stronger clustering of 
spin vectors around the Cartesian mesh basis vectors than 
for the high-resolution run. Spurious grid alignment is thus 
present at z — in both simulations to some degree. 

A4 Alignments with the Tidal Field 

In Figure [A"2l we show the scale dependence of the alignment 
between the galaxy stellar, gas and dark matter angular mo- 
mentum and the three eigenvectors of the tidal field for the 
low-resolution run. These results should be compared with 
results for the high-resolution simulation given in Figure [TT1 
We find that the results are perfectly consistent for all sam- 
ples except the low-density low-mass bin where we find a 
slightly stronger anti-alignment with the intermediate prin- 
cipal axis of the tidal field in the low-resolution run. Also, 
the weak anti-alignment with the third principal axis is not 
seen in the high-resolution run. It is likely that this differ- 
ence for low-mass galaxies in low density regions is due to the 
anisotropy induced by the Cartesian mesh which is expected 
to be stronger for the low-resolution run. We conclude that 
the resolution of our high resolution run is still insufficient 
for galaxy orientations to be not affected by spurious grid 
relaxation in low density regions at low redshifts. It is thus 
not clear whether the alignment of spins with the interme- 
diate axis that we find at higher redshifts would persist in 



low-density environments also down to z — in the absence 
of numerical effects. This question has to be addressed with 
even higher resolution simu lations than th e ones we used in 
our analysis. The results of iNavarro et alj ((2004), who find 
a survival of alignment with the intermediate axis for the 
four isolated galaxies they simulated, is however indicative 
that such alignment should persist to some degree. 
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Figure A2. Same as Figure[9]but for the simulation with two-times lower spatial resolution at redshift z = 0. The corresponding results 
for the high-resolution run are shown in Figure [TT1 



